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L Introduction 


r In  this  paper  an  analysis  is  made  of  the  reflection  of  a pulse  oi  sound  from 
a ooint  source  in  a non-abeorbing  liquid  plane  boundary,  under  the  circumstance 
that  the  index  of  refraction  is  <1.  The  analysis  is  equally  applicable  to  the 
total  reflection  of  electromagnetic  radiation  from  a vertical  dipole,  in  the  absence 
of  absorption.  The  results  clarify  the  mechanism  of  propagation  of  the  head 
wave  which  has  been  observed  in  acoustic  experiments.  It  is  shown  without 
approximation  that  % pulse  is  received  in  the  first  medium  at  a time  correspond- 
ing to  propagation  along  the  surface  with  wave  velocity  of  the  second  medium. 

Under  certain  circumstances  a pole  in  the  complex  plane  of  integration 
maW-w  a contribution  to  the  field  at  long  distances.  However  the  pole  always 
lies  on  the  imaginary  axis  for  non-abeorbing  media,  giving  rise  to  an  exponen- 
tially decreasing  term  at  long  distances.  Thus  no  true  (undamped)  surface 
waves  can  occur.  The  contribution  from  this  pole  is  absent  in  the  electro- 
magnetic case  when,  as  is  usual,  the  magnetic  permeabilities  of  the  two  media 
are  equal 


Figure  I shows  a point  source  at  P,  a height  h above  an  infinite  plane 


'Contribution  from  the  Scrippe  Institute  of  Oceanography.  This  work  represents  one  of 
the  results  of  research  earned  out  under  coo  tract  with  the  Bureau  of  Ships,  Navy  Department 
t Present  address:  Department  of  Phyrics,  University  of  Pittsburgh. 
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surface  separating  two  non -absorbing  media.  The  wave  velocity  is  greater  in 
medium  2 than  in  medium  1.  The  signal  is  received  at  S,  a height  t above  the 
bounding  surface. 

Let  the  source  be  harmonic,  of  circular  frequency  u.  In  the  acoustic  case, 
when  both  media  are  liquids,  the  acoustic  potentials  4,  and  satisfy  in  their 
respective  media 


(1) 


A$,  — 0, 

A*,  + ***•  - 0 


with  the  boundary  conditions,  on  the  plane  t — 0 


(2) 


& 

dx  dt‘ 


In  these  equations  p denotes  density,  e wave  velocity,  and  k — »/ c. 

In  the  electromagnetic  case,  when  both  media  are  non-conducting,  and  the 
source  is  a vertical  dipole,  the  Herts  potentials  t,  end  r,  satisfy  [1] 


(3) 


A»,  + kWi  - 0, 
Ati  + ■*  o 


with  the  boundary  conditions,  on  the  plane  * - 0 


(4) 


dr,  9rt 
dz  ~ dt  • 


In  equation  (4)  p is  the  magnetic  permeability.  It  will  be  noticed  that  when 
Hi  — p,  the  electromagnetic  problem  is  less  general  than  the  acoustic,  since  the 
ratio  of  the  coefficients  of  t,  and  rt  in  the  first  of  equations  (4)  is  k\/k\  — cj/c? , 
whereas  p,/p,  is  independent  of  c,/ci . 

The  field  in  medium  2 has  been  discussed  previously  [2],  in  this  paper  the 
field  in  medium  1 will  be  analysed.  In  the  main  our  results  confirm  those  of 
previous  investigators  [3]  but  some  new  features  are  introduced,  to  wit:  (1) 
The  time  of  arrival  of  a pulse  is  evaluated  exactly,  and  shown  to  agree  with  the 
conclusions  inferred  from  the  approximate  evaluation  of  the  field  from  a harmonic 
source,  (t)  The  validity  of  the  method  of  steepest  descents  in  this  problem  is 
examined  by  evaluating  the  next  higher  approximation,  and  a number  of  the 
conditions  for  the  validity  of  the  method  are  shown  to  have  simple  and  natural 
physical  interpretations.  The  amplitude  dependence  of  the  head  wave  is  also 
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given  a ample  geometrical  interpretation.  The  head  wave  denotes  the  wave 
which  in  the  region  of  total  reflection  arrives  before  the  reflected  wave,  and  may 
arrive  before  the  directly  transmitted  wave.1  (5)  Under  certain  circumstances 
a pole  lying  on  the  imaginary  axis  in  the  complex  plane  of  integration  makes  a 
contribution  to  the  field,  yielding  an  exponentially  damped  wave  at  long  dis- 
tances in  medium  1.  The  circumstances  are  such  that  the  pole  makes  no  contri- 
bution in  the  electromagnetic  case  when  m,  *=  The  contribution  of  this  pole 
does  not  seem  to  have  been  previously  remarked.* 

II.  Notation  and  Basic  Formulas 

We  employ  our  previous  notation:  the  acoustic  potential  at  5 (Figure  1) 
with  the  coordinates  (x,y,z)  is 

(5)  *.-*,  + *. 

R is  the  distance  PS  (Figure  1),  and  the  reflected  wave  4>r  is 

(6)  4>r  « ik\  f du  - i f,Ja(ktru)  exp  [tfc,(*  + h)Vl  - u*|, 

Jo  V 1 — u 

rn  t flVl  - tt*  — V a — u* 

l)  /r  " fi  vT^P  + 

In  equations  (6)  and  (7),  r — (x*  + a — c,/cj  is  assumed  <1,  0 equals 

Pi/p,  (acoustic)  or  p, a/n,  (electromagnetic),  and  may  have  any  positive  value. 
Equation  (7)  for  the  reflection  coefficient  /,  is  equation  (10)  of  reference  2. 
Equation  (6)  follows  readily  from  equation  (6)  of  reference  2 by  the  procedure 
used  to  obtain  equation  (13)  of  reference  2.  Details  of  the  derivation  of  equa- 
tion (6)  are  also  given  by  Ott  |3],  The  correct  behavior  of  the  solution  at  in- 
finity  in  either  medium  is  guaranteed  by  choosing  Vo- t?  and  y/\  — u* 


‘Such  early  arriving  waves  have  been  observed  in  geophysical  prospecting  by  explosive 
sound*  (cf.  Myokat,  ref.  3),  and,  in  underwater  propagation  by  J.  L.  Worxel  and  M.  Ewing, 
EzjpUmon  sound*  in  thallow  water,  Geol.  8oc.  of  Am.,  Memoir  27,  1948. 

•Of  the  authors  listed  in  reference  3,  Lamb  solved  not  the  problem  under  consideration 
• here,  but  some  related  problems.  Jeffreys  uses  operational  methods  throughout.  Muskat 
obtained  a valid  formal  solution,  but  evaluated  the  resultant  integrals  by  means  of  a relatively 
crude  approximation.  Ott’s  analysis  parallels  our  own  but  be  does  not  follow  the  steepest 
deeoeot  contours  in  the  complex  plane  with  the  care  which  we  have  used.  In  essence  our 
evaluation  of  the  integrals  oon firms  Ott’s  results,  but  extends  them  by  including  the  small  but 
theoretically  interesting  contributions  from  the  pole.  Kruger  evaluated  the  integrals  without 
uring  the  method  of  steepest  descents,  for  the  electromagnetic  case.  However,  his  analysis 
is  immediately  applicable  to  total  reflection  only  when  k - 0,  and  uses  the  restrictive  assump- 
tion that  mi  — mi  • Brehovtkih’s  work,  in  Russian,  seems  to  parallel  our*,  but  he  makes  no 
mentioo  of  the  pc'.e. 
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positive  real  for  small  positive  u and  positive  imaginary  for  large  positive  u. 
The  choice  of  sign  is  achieved  by  drawing  the  cuts  through  a and  1 upward,  as 
shown  in  Figure  2.  It  will  prove  convenient  to  draw  the  cuts  through  —a  and 
— 1 downward  as  shown;  for  the  purposes  of  this  section  their  directions  arc 
not  significant.  The  contour  of  equation  (6)  runs  from  0 to  « on  the  real  axis 
without  crossing  any  of  the  cuts.  For  z+h  > 0 

exp  \ikxy/r*  + (z  4- 

Vra  + (z  + A)* 

(8) 

- f dt  -o (ri)  exp  {-(*  + h) 

Jo  vi  - k : 


Equation  (8)  is  well  known;  it  implies  that  the  real  part  of  Vr— ~k*,  is  always 
>0  on  the  contour,  (cf.  8tratton  (bj)  or  equivalently,  that  the  contour  passes 
below  the  singularity  at  t — k,  (cf.  Watson  [6]).  Letting  t — kxu  in  equation 
(8),  remembering  ifc,  > 0,  and  noting  that  vV  -1-  1 * — * V 1 — u*  is  con- 
sistent with  the  definitions  of  the  signs  of  the  radicals  in  equations  (6)  and  (8), 
it  con  be  seen  that 


(9)  i 


(10) 


_ eik'4  [~ fl(*  + h)/d  - V^?l 

+ h)/d  + Va  - r*J 

= 2ikt  f du  ~A~ i F(M)Jo(kx ru)  exp  {fk,(z  + h)V  1 - u*), 

Jo  VI  - u 


(11) 
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(12)  F(u) 


iSVl  - r*  + Va  - r*  0 y/ 1 — ua  + y/cf — u*’ 


(13)  d = vV  + (2  + JO*  , r-r/d,  Vl  - t * » 


The  distance  between  the  receiver  and  the  image  of  the  source  is  d = P’S 
(Figure  1);  r is  the  3ine  of  the  angle  of  incidence  TPP'  (Figure  1).  In  equation 
(7)  u equals  the  sine  of  the  angle  of  incidence  (see  [2]);  equations  (6)  and  (11) 
are  integrals  over  real  and  complex  angles  of  incidence.  The  amplitude  and 
phase  of  <£.  are  precisely  those  expected  for  z spherical  wave  propagating  along 
the  (geometrically  reflected  ray  PTS  (Figure  1),  if  we  assume  this  wave  to  be 
reflected  with  the  plane  wave  coefficient  of  equation  (7),  with  u >=  r.  The  sign 
of  y/ a — r 1 is  identical  with  that  of  y/a  — u1  when  u = r.  We  infer  therefore 
that  all  deviations  from  ray  or  geometrical  propagation  are  contained*  in  the 
term  4,  . 

The  poles  of  the  integrand  of  equation  (11)  lie  at  the  roots  of 


(14) 
or  at 


(15) 


0y/l  - u*  + Va  - u = 0 


1 <0, 


a < 0 < 1, 


0 < or. 


In  that  part  of  the  complex  plane  which  is  above  the  real  axis,  we  denote 
by  I the  region  to  the  left  of  the  cut  at  a,  by  II  the  region  between  the  cuts  at 
a and  1,  and  by  III  the  region  to  the  right  cf  the  cut  at  1.  Since  the  cuts  may 
(and  will  later)  be  drawn  so  that  they  intersect  the  positive  imaginary  axis,  we 
further  distinguish  between  region  I,  , containing  thoee  points  in  region  I which 
lie  in  the  first  quadrant  of  the  complex  plane,  and  I,  , containing  those  points 
in  region  I which  lie  in  the  second  quadrant,  and  similarly  for  regions  II  and 
III.  The  fourth  quadrant  is  designated  by  IV.  Some  of  these  regions  are 
shown  in  Figure  2.  In  regions  I, , I, , II,  , III,  , and  IV  the  real  and  imaginary 
parts  of  (1  — u*)'",  (a*  — u*)‘  *>  and  u/(l  — u*),/*  are  negative  in  the  following 

This  formulation,  in  which  — 0 in  the  geometrical  limit  Jfc,  — ► «=  appears  more  logical 
than  the  customary  procedure  of  ignoring  the  dependence  on  a,  0,  and  angle  of  incidence  of  the 
reflection  coefficient,  whereby  the  field  in  medium  1 is  regarded  as  the  sum  of  spherical  waves 
radiating  from  sources  of  equal  strength  at  P and  P'  (Figure  1),  phis  s contour  integral,  which 
integral  is  not  sero  in  the  geometrical  limit.  Cf.  Stratton,  op.  oil.,  p.  573  ff.,  or  Ott,  op.  ext. 


78 


EDWARD  OERJUOY 


regions: 


(Re  (1  - 

- u3y/3 

< 

0: 

IN, 

6m  (1  - 

- uyn 

< 

0: 

I,  , 

II, 

(Re  (a*  - 

- u’)'" 

< 

0: 

II, 

N;  , III, 

6m  (a*  - 

- uy/2 

< 

0: 

I,  , 

II, 

<Reu/(l  - 

- uy' * 

< 

0: 

I*. 

II,  , IN, 

6m  u/(l  - 

- uyn 

< 

0: 

Ill, 

,iv. 

Otherwise  they  are  positive.  The  real  and  imaginary  parts  of  the  radicals 
change  sign  discon tinuously  at  the  cuts  but  become  zero  on  the  real  or  on  the  im- 
aginary axis  if  they  change  sign  crossing  it. 

Equations  (16)  show  there  are  no  roots  of  equation  (14)  in  region  IV.  Con- 
sequently, as  shown  for  similar  integrals  by  Muskat  and  Ott  [3]  we  obtain 

(17)  4>,  = f du  U— F(u)Ho\kiru)  exp  |tfc,(z  + h)Vl  - u1). 

Jc  vl  — u 

The  contour  C of  equation  (17)  runs  from  i =>  to  «=  on  the  real  axis,  as  shown  in 
Figure  2.  Replacing  Ho  ’ by  its  asymptotic  expansion,  equation  (17)  becomes 

(18)  *,  - (^’V"  Jc  du(r-^)'',F(u)  exp 

(19)  A{u)  = (t  + h)  Vl  — a*  + ru. 

In  equation  (17)  — x < arg  u < x.  To  prevent  circling  the  singularity  at  the 
origin,  the  cut  at  —a  is  extended  to  the  origin  along  the  negative  real  axis  as 
indicated  by  the  stippling  in  Figure  2. 

We  evaluate  4>,  , equation  (18),  by  the  method  of  steepest  descents.*  There 
are  two  distinct  cases:  r < a,  considered  in  the  immediately  following  section 
and  t > a,  which  is  deferred  to  section  IV. 

III.  Solution  for  x < a 

The  saddle  points  of  the  integral  (18)  lie  at  the  roots  u = ±r  of  the  equa- 
tion A’(u)  •»  0,  with 

(20)  A’(u)  = r - %±ML. 

vl  — u 

We  shall  be  concerned  only  with  the  point  u — /,  which  for  t < a lies  in  I,  . 
A(u)  is  pure  real  on  the  real  axis  between  0 and  1,  and,  in  It  , has  a maximum 
at  v **  r.  Along  the  contour  (Re  A(u)  *■  A(r)  through  u — r in  I,  , 6m  A(u) 


•Wataon,  op.  cit.,  p.  2. 
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can  have  no  maxima  and  but  a single  minimum,  at  u — r.  Evaluation  of  4>, , 
equation  (18),  by  the  method  of  steepest  descents  now  proceeds  in  the  usual  way. 
To  completely  justify  the  analysis  however  it  is  necessary*  to  make  sure  that  it 
is  possible  to  deform  the  contour  C of  Figure  2 into  the  contour  <Ke  i4(u)  — X(r) 
passing  through  u «*  r in  I,  . 

The  argument  required  is  similar  to  that  employed  previously.  Use  is 
made  of  equations  (16)  and  (20),  together  with  the  remark  that  the  contour 
is  parallel  to  the  imaginary  (real)  axis  if  and  only  if  A'(u)  is  pure  real  (imaginary). 
If  u — pe'*,  then  in  III,  and  IV,  for  large  p 

(21)  /4(u)  =*  p{[r  cos  0 (z  + h)  sin  0]  + t[(z  -j-  h)  cos  0 + r 6in  0] } 
while  in  I,  and  I, , fcr  large  p 

(22)  /4(u)  = p{  [r  cos  0 + (z  -j-  h)  sin  0]  + t'[r  sin  0 — (z  4-  h)  cos  0]} . 

It  can  be  seen  from  equation  (16)  that  the  only  roots  of  equation  (14)  which  can 
possibly  concern  us  lie  in  region  II.  We  can  conclude  that  the  contour  C can 
be  deformed  into  the  contour  <R e A (u)  — A(r),  which  we  denote  by  C’  (Figure  3), 
without  crossing  any  poles. 

We  recall  that  u is  to  be  interpreted  as  the  sine  of  the  angle  of  incidence  of 
an  arbitrary  ray  from  the  source.  The  equation  u — r for  the  saddle  point  of 


Fioure  3 

equation  (18)  corresponds  therefore  to  the  assertion  that  the  deviations  from 
ray  propagation  arc  in  this  case,  r < a,  determined  primarily  by  radiation  along 
the  direction  of  the  expected  geometrical  path.  This  interpretation  is  confirmed 
by  the  fact  that  exp  jtA,rl(T)}  * exp  is  a phase  factor  common  to  <t>.  , 

equation  (10),  and  to  this  estimate  of  4>t  • Moreover  since  F{u)  = 0 at  u = r, 

This  point  has  not  always  been  sufficiently  emphasised. 
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4,  is  seen  to  t-e  small  compared  to  4,  (in  the  limit  of  high  frequencies  and  large 
distances,  where  the  method  of  steepest  descents  is  applicable).  Theee  results 
are  in  complete  agreement  with  those  of  Ott  (3). 

IV.  Solution  for  x > a 

When  r > a it  is  necessary  to  deform  the  contour  into  the  shape  C,  followed 
by  C,  of  Figure  4 in  order  to  obtain  a convergent  contour  integral  through  the 
Baddle  point  u =*  r.  We  write 

(23)  4.  = 4,(Ci)  + 4,(C,) 

where  4,  is  given  by  equation  (17),  and  both  4, (C,)  and  4,(6’,)  are  given  by 
equation  (17)  with  the  substitution  of  C,  or  C,  for  C.  Replacing  Hol)  by  its 
asymptotic  expansion,  4,(C,)  is  given  by  equation  (18)  with  C,  instead  of  C. 


Since  Va1  - u*  does  not  occur  in  4(u)  the  argument  in  the  preceding  section 
concerning  the  determination  of  the  contour  (Re  A (u)  = A(t)  is  wholly  un- 
changed, except  for  the  fact  that  the  contour  extends  to  infinity  on  the  left  in 
II,  instead  of  I, . The  only  question  to  be  settled  is  whether  deformation  of  the 
contour  requires  crossing  any  of  the  poles  in  II.  This  question  will  be  considered 
below.  Except  for  possible  contributions  from  poles  in  II  the  integral  4,(C,) 
is  seen  to  yield  precisely  the  same  results  as  did  4,  in  the  preceding  action.  The 
contour  into  which  C,  is  deformed  is  termed  CJ . 

We  proceed  to  4,(C1):  The  first  term  in  F(u),  equation  (12),  is  merely  a 
constant.  Referring  to  equation  (17)  and  to  Figure  4,  it  is  evident  that  in  4,(C,) 
the  contour  C,  can  be  closed  at  infinity  for  this  first  term  in  F(u),  since  the  integral 
containing  this  term  contains  no  terms  involving  \/a*  — u*.  Nor  does  it  contain 
any  poles.  Consequently  the  first  term  in  F(u)  results  in  a vanishing  integral 
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I 


-1 f 

J c, 


u 


du 


Va2  - u2 


Vl  - u2  0 Vl  - u2  + Vo5  - u2 


and  we  have 

*.(C0 

(24)  

• exp  (t*,(z  + h)VT-  u2)Hol){kiru). 

In'll,  the  contour  C,  lies  along  the  cut.  The  contour  C,  can  be  deformed  so  a a 
to  lie  along  the  cut  in  I, , since  there  are  no  poles  in  I, . Because  the  integrand  in 
equation  (24)  vanishes  at,  u — a the  integral  around  a small  circle  about  u — a 
is  zero.  It  follows  that 

A II  I H _ */•!•  ■ /1*l  ■ 

(26) 


♦.(CO 


-2,7:,  / 


du 


O.  u“  **(l  - u2)  - (^ 


• exp  |t*,(z  -4*  h)  Vl  - u*{ Hi1 ’(fc.ru). 

In  equation  (25)  the  contour  C»  runs  in  IIt  along  the  cut,  starting  at  u — a 
and  going  to  infinity. 

We  desire  the  “best”  contour  for  the  branch  line  integral  $,(C,);  namely 
a contour  which  appears  likely  to  minimize  the  error.  Although  no  saddle 
point  can  be  found  for  $,(C,),  it  is  very  reasonable  to  deform  C,  into  a contour 
along  which  the  phase  of  the  integrand  remains  constant.  Such  a steepest 
descent  contour,  whether  or  not  it  passes  through  a saddle  point,  still  seems 
least  likely  to  introduce  complicating  cancellations  and  reinforcements,  par- 
ticularly in  the  large  k,r  limit  (in  which  we  are  ultimately  going  to  be  most 
interested). 

Consequently  an  approximate  expression  for  $,(C,)  is  obtained  much  as 
before  by  replacing  J/J1’  in  equation  (25)  by  its  asymptotic  expansion,  and 
determining  the  contour  <R*  A(u)  — A(^)  into  which  the  cut  through  a must  be 
deformed.  Only  if  this  deformation  of  the  contour  requires  crossing  any  of  the 
poles  in  II  must  their  contributions  be  considered.*  The  contour  (Re  A (t<)  ■»  A (a), 
which  we  term  C[ , is  shown  in  Figure  5 where  the  cut  through  a has  been  de- 
formed so  that  it  coincides  with  C|  . The  deformation  is  justified  as  in  the 
preceding  section.  The  integral  along  C(  i8  computed,  as  in  the  method  of 
steepest  descents,  by  expanding  the  integrand  about  the  point  u — a.  The 
final  result  is  an  expression  for  the  head  wave: 

2ai 

(26)  *,(C{)  *** 


*,0(1  -_^T7i  r,/2(r(l  - a*) 1/1  - (z  + h)a\ 


i/i 


and 

(27) 


* - {z  + h)  Vl  - a1  + 


m. 


*It  is  definitely  not  correct  to  try  to  evaluate  the  integral  along  an  arbitrary  contour 
between  a and  <*>,  The  exact  value  of  the  integral  is  independent  of  the  choioe  of  contour,  of 
oouree.  But  because  the  integral  is  not  evaluated  exactly,  the  estimated  value  of  the  integral 
can  depend  on  the  choice  of  contour,  as  apparently  ’ appena  in  the  Sommerfeld  problem. 
T.  Kahan  and  Q.  Eckert,  Jour.  Phya.  Rad.,  10,  165  (1949)  and  Phys.  Rev.,  76,  406  (1949). 
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Equation  (26)  shows  that  in  the  limit  of  high  frequencies  and  large  distances  the 
amplitude  of  $,{£[)  is  small  compared  to  4>.  , equation  (10).  This  was  to  be 
expected,  since  4>,(C[)  is  part  of  the  deviation  from  geometrical  propagation. 

We  return  to  the  problem  of  whether  or  not  poles  are  crossed.  It  is  apparent 
from  the  above  discussion  that  we  might  just  as  well  have  deformed  the  contours 
to  C[  and  C't  before  replacing  Hlau  by  its  asymptotic  expansion.  Figure  5 
illustrates  some  intermediate  stages  in  the  deformation  of  the  contour  C of 
Figure  2 into  integrals  along  C[  and  CJ  and  shows  the  cuts  through  u — —a  and 
u «*  1.  The  contour  C[  coincides  with  the  cut  through  a and  is  drawn  stippled, 
as  is  the  cut  from  u — — a to  u — 0.  The  portion  of  C't  to  the  left  of  u — r is 
also  shown.  Pegions  I,  and  I,  lie  below  and  to  the  left  of  C{ . In  Figure  5 the 
contour  C has  been  bent  around  into  region  II  between  C[  and  Ci  , but  has  not 
yet  been  extended  to  infinity.  The  dotted  line  is  a deformation  of  C which 
has  not  crossed  the  imaginary  axis. 


Referring  to  equation  (15)  it  is  clear  that  if  1.  < 0 it  is  possible  to  extend 
to  C in  Figure  2 to  the  dotted  C in  Figure  5 to  the  solid  C in  Figure  5 and  so  on 
to  infinity  in  II, , without  crossing  any  poles.  When  (i  < a the  poles  lie  on  the 
real  axis  between  —a  and  a in  region  II.  Thus  C can  again  be  deformed  to 
infinity  in  II  without  crossing  a pole.  For  a < 1 3 < 1 the  poles  lie  on  the  im- 
aginary axis;  the  only  root  of  significance  is  the  positive  imaginary  root  in  II. 
Let  u — iyx  (ty,)  at  the  intersection  of  C[  (CJ)  with  the  imaginary  axis.  If  the 
pole  lies  at  y , > y,  (Figure  5),  the  contour  can  obviously  be  deformed  from  the 
dotted  C to  the  solid  C without  crossing  the  pole,  and  so  on  to  infinity.  If  the 
ordinate  of  the  pole  is  <y,  , the  pole  lies  below  the  cut  C{  , and  is  again  never 
crossed.  The  only  case  that  remains  is  a < 0 < 1,  y,  < y,  < y» , as  illustrated 
in  Figure  5.  In  this  case  the  contour  cannot  be  extended  from  the  dotted  C to 
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the  solid  C without  crossing  the  pole,  so  that  to  the  integral  along  the  solid  C 
must  be  added  the  integral  around  the  closed  loop  surrounding  y,  ; this  integral 
readily  yields 

*.(P)  - 4 ik,f  exp  {*,[(*  + h)(l  - aTVU  - 0*)!/,]| 

(28)  (1  ~ 

$,(p)  is  the  contribution  from  the  pole  and  in  this  case  where  a < $ < 1, 
Vi  < V,  < Vi  we  have 

(29)  ♦,«?,)  - *,(p)  + ♦,(C0 

where $,(C()  is  still  given  by  equation  (20).  In  equation  (28),  K„  is  the  modified 
Hankel  function.  The  condition  yt  < y,  < yt  can  be  written  as 

(a®  (i  - »■)•'■  + < {rz$T  < rfv 

We  have  been  unable  to  interpret  $,(p)  physically.  Since  for  large  z 


4>,(p)  decreases  exponentially  for  large  r,  and  therefore  can  hardly  be  interpreted 
as  a surface  wave  in  this  problem  which  neglects  absorption,  despite  the  r_,/* 
spreading  factor  inferred  from  equations  (28)  and  (31).  It  is  likely  that  at  large 
distances  the  magnitude  of  4,(p)  is  small  compared  to  the  error  made  in  deriving 
equation  (26).  It  is  possible  that  $,(p)  is  a portion  of  the  difference  between  the 
accurate  and  asymptotic  value  of  $,(C,)f  or  it  may  be  a spurious  result  of  the 
analysis.  All  that  can  be  said  is  that  our  choice  of  "best”  contour  C{  leads 
naturally  to4,(p)  when  a < 0 < 1 and  equation  (30)  is  satisfied.  In  the  electro- 
magnetic case  when  p,  « m*  , we  have  8 = a,  so  that  a < 8 < 1 is  not  possible 
under  these  circumstances. 

V.  Geometrical  Interpretation 

In  this  section  we  assume  the  head  wave  propagates  from  P to  S along  the 
path  PABS  (Figure  1)  with  velocity  e,  Along  AB,  and  show  that  this  assumption 
leads  to  simple  geometrical  interpretations  of  the  results  of  the  preceding  sections. 
In  Figure  1 

(32)  AB  = r - + ^ . 

Vl  - a 

Musks t [3]  has  pointed  out  that  the  phase  difference  between  P and  S corresponds 
to  the  assumed  propagation,  i.e.,  referring  to  Figure  1,  that 

(33)  kt(PA)  + aki(AB)  + 1 k,(BS)  = 


i 
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holds  with  i given  by  equation  (27).  He  also  remarks  that  PABS  is  a path  of 
least  time,  in  the  sense  of  Fermat’s  principle.  The  surfaces  of  constant  phase 
— constant  intersect  the  1,2-plane  of  Figure  1 in  lines  parallel  to  USW  and 
perpendicular  to  BS.  The  wave  fronts,  obtained  by  rotating  USW  about  the 
s-axis,  are  obviously  cones  with  axis  OP.  These  cones  have  an  angle  «•  CSW  — 
the  complement  of  the  critical  angle. 

Consider  the  two  neighboring  rays  BS  and  B'S'  intersecting  the  wave 
front  USW  (Figure  6).  Rotate  the  section  BSS'B'  about  OP.  All  the  energy 


lying  within  the  volume  formed  thereby  which  in  unit  time  crosses  the  surface 
formed  by  rotating  SS ' ultimately  crosses  the  surface  formed  by  rotating  VV’, 
V and  V'  being  intersections  of  BS  and  B'S'  with  the  neighboring  plane  of 
constant  phase  U'VW’.  Calling  the  pressure  at  S p(r,z)  and  at  V p(r+drj+dt:) 
and  noting  that  SS'  = VV'  we  have 

(34)  2wt p*(r,  2)  = 2r(r  + dr)p*(r  + dxjt  + dt). 

Equation  (34)  implies  that  for  points  on  the  same  ray  BS,  p is  proportional  to 
r_,/a  and  independent  of  2.  This  is  precisely  the  result  yielded  by  equation 
(26)  from  which,  using  equation  (32),  $4(C 0 is  proportional  to  (AB)~*/*r_,/\ 
■ Clearly  equation  (26)  becomes  grossly  inaccurate7  as  AB  approaches  zero 

I length,  since  the  solution  must  be  finite  at  all  points  other  then  P.  In  Figure 

j 6,  S"  is  a point  in  the  field  receiving  the  head  wave  via  the  ray  B"S".  B"M  is 

the  wave  front  starting  at  B".  It  terminates  abruptly  at  M,  since  the  critical 
ray  PA  is  reflected  into  AML  and  since  there  is  no  refracted  wave  along  the 
surface  to  contribute  to  the  head  wave  at  points  on  the  surface  cloeer  to  the 
origin  0 than  A.  We  may  expect  our  assumption  of  propagation  along  PABS 
to  be  valid,  and  equation  (26)  to  be  a good  approximation,  as  long  as  B"M 
subtends  a large  number  of  Fresnel  zones  constructed  from  point  S".  In  other 

’Brehovskih,  see  (3],  by  a different  method  obtained  instead  of  equation  (26)  a modified 
expression  which  he  claims  accurately  represents  4,(C{)  for  email  or  sero  values  of  AB. 


I 

1 

! 


\ 
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words  it  is  necessary  that 
(35) 


S"M  - S"B" 


» 1 


when  S"B"  » MB")  equation  (3-5)  simplifies  to 
(36)  (1  - a*)(AB)'  » 2Xr/Vl  - a*. 


Equation  (26)  is  the  loading  term  in  an  asymptotic  expansion  of  6,(C,)t 
presumably  valid  in  the  limit  as  kxr  — *®.  For  finite  kxr  the  ratio  of  the  next 
term  in  the  asymptotic  expansion  to  the  leading  term,  equation  (26),  furnishes 
an  estimate  of  the  error  made  in  using  equa  ion  (26)  to  represent  the  head  wave. 

aa*  yMuy 

(37)  «'•(*, rv)  - (^)‘"  exp  (.('  .ru  - x/4))(l  + 


where  | A | < 1.  If  in  equation  (25),  we  ignore  the  term  in  A,  introduce  the  new 
variable  17  which  runs  from  0 to  00  as  u traverses  C[  from  a to  <*> , and  expand 
the  terms  in  the  integrand  exclusive  of  the  exponential  terms  in  a power  series, 
we  obtain 


(38) 


*,(CQ 


4iak\n<?'' 

- a)1'*  rl/M  1 - a')l/1  - (z  + h)a]i/J 

• f d,e-'*‘V/,(l  + **«  + •••). 

Jo 


The  leading  term  in  equation  (38)  yields  equation  (26).  In  equation  (38) 
i — X),  bt , j -»  1 to  4,  where  the  St  are  defined  as  follows:  For  j — 1 to  3, 

(39)  Qf  = BJ!  + «,(*>)  + •••] 

the  functions  Qt  , j — 1 to  3,  being  respectively  u,/*,  [/3*(1—  w*)  — (a*— u*))_l, 
and  —idu/dij.  i«  is  defined  by 

(40)  (a*  - uT’  « Bj'W  + a«(tij)  +.-•]. 

In  equations  (39)  and  (40)  the  B( , j °*  1 to  4,  are  numerical  coefficients. 

In  equation  (38)  the  absolute  value  of  the  ratio  of  the  term  in  6 to  the  leading 
term  is  Zb/(2kxd),  which  must  be  «1  if  equation  (26)  is  to  be  accurate.  The 
expression  Z6/(2kxd)  where  two  or  more  of  the  ZSt/(2kxd)  are  large  can  sometimes 
be  made  «1  by  cancellation;  such  cancellation  is  seldom  possible  in  the  omitted 
term  of  equation  (38),  of  order  ij*  in  the  integrand,  whose  ratio  to  the  leading 
term  will  be  ~5 *(&,<£)*.  We  conclude  therefore  that  it  is  probably  necessary 
that  each  of  the  Zh{/{2kxd)  be  «1,  in  order  that  equation  (26)  be  accurate. 

The  inequality  corresponding  to  j = 3 is 

3(r  -f-  k) . 

2(1  - a*),/Ji fc,(A£)J  <<[  ** 


•Wstaon,  op.  at.,  p.  219. 
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It  is  a confirmation  of  our  geometrical  interpretation  that  equation  (41)  ia 
equivalent  to  equation  (36).  Apparently  the  geometrical  approximationa  in- 
volved in  deriving  equation  (36)  implied  z » h.  That  the  right  hand  aide  of 
equation  (36)  ahould  depend  on  (z  -f  h)  can  be  inferred  from  equationa  (6)  and 
(11)  which  ahow  that  4>r  and  4>,  are  function  of  z and  h only  through  the  sum 
(z  + h). 

We  may  add  that  in  equation  (25)  the  contribution  of  the  term  in  A can 
be  estimated  by  expending  that  term  in  a power  series  about  u — a,  leading  to 
the  inequality  8 kxra  »1,  a typical  condition  for  the  validity  of  geometrical 
propagation.  Some  of  the  inequalities  3i,/(2fc,d)  «1  imply 


(42) 


4 aki[AB) 


l « 1. 


We  have  not  interpreted  equation  (42).  It  can  be  inferred  from  35</(2l:,of)  « ! 
that  (1  — a*)  and  $ must  not  be  too  close  to  zero,  conditions  to  be  expected 
from  equation  (26),  since  0,(C,)  cannot  increase  indefinitely  as  (1—  a1)  and  0 
approach  zero. 


VI.  Solution  for  Complex  to 

To  this  point  we  have  been  concerned  with  positive  real  kx  and  k2  . To 
determine  the  received  signal  resulting  from  a pulse,  we  require  the  solutions  for 
fomnlov  vaIiip*  of  <.i  where  8m  oi  > 0.  We  retain  kx  — w/c,  , k2  — w/c,  , cx  and 
c,  real.  The  previously  adopted  means  of  obtaining  the  solution,  based  on  the 
expansion  of  e'kM/R  in  plane  waves  as  first  used  by  Weyl  [7],  would  in  the  present 
case*  involve  divergent  integrals.  We  can  follow  Sommerfeld  [8]  however; 
by  analytic  continuation  it  can  be  seen  that  equation  (8)  is  correct  with 
(A  + z)  > 0,  in  the  extended  range  r > 0,  8m  ki  > 0,  where  (Re  y/  ? — k]  > 0. 
The  value  of  y/ 11  — k\  at  points  ? <k\  ,kx  real,  is  determined  by  continuation 
of  y/ 1*  — k\  from  the  upper  half  plane  0m  kx  > 0.  We  note  that  (Re  y/ 1*  — k* 
cannot  change  sign,  as  fc,  is  varied  continuously  for  fixed  l,  unless  8m  kx  becomes 
zero,  and  that  the  integral  therefore  converges,  with  Vf*  — k*  as  defined,  for 
r *=  0 and/or  k , — 0. 

Using  equation  (8),  with  ( h — z)  replacing  (A  + z),  so  as  to  represent  , 
equation  (5),  we  can,  as  does  Sommerfeld,  obtain  a solution  for  complex  fc,  . In 
fact,  letting  t — kxu  in  equation  (8)  and  the  other  relevant  integrals,  the  solution 
in  medium  I is  given  by  equations  (5)  and  (9)— (13),  except  that  the  contour  for 
4,  , equation  (II),  now  runs  from  0 to  ® along  the  line  arg  u — — arg  A,  . The 
cuts  through  u = 1 and  u — a,  formerly  through  t = kx  and  ak,  , approach 
infinity  along  arg  v *»  (x/2)  — arg  fc, . The  cuts  through  u — — 1 and  u — — a, 

The  plane  wave  expansion  of  eikM/R  appears  to  us  to  involve  divergent  integrals  in 
medium  2 whenever  kt  is  complex  even  with  real  Jt,  . When  k,  is  complex,  for  any  x and  y, 
there  are  values  of  the  azimuth  angle  <!>  in  equation  (12)  of  reference  2 for  which  the  integrand 
becomes  infinite  as  9 approaches  its  upper  limit  of  t/2  — * «=.  This  difficulty  does  not  seem  to 
have  been  remarked  by  later  writers,  e.g.  Ott 
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which  were  downward  in  the  f-plane,  approach  infinity  alo.ig  arg  u — — (v/2) 
— arg  k{ . With  these  provisos,  equation  (11)  can  be  seen  to  represent  an  analytic 
function  of  fc,  in  the  entire  range  dm  k , > 0;  as  always,  h,  — afc,  , a < 1.  This 
result  is  most  important.  Knowing  the  solutions  are  analytic  functions  of  kx 
enables  us,  as  in  the  next  section,  to  extend  integrals  over  u into  the  complex 
u>-plane.  If  kx  =•  pc’*,  0 < 8 < r,  the  behavior  of  <t>,  , equation^ll),  for  fixed  r 
and  2 as  p — »<*>  is  determined  by  exp  | — p(/t-f 2)[sin  8 (Re  y/l  — u*  + cos  8 
dm  y/l  — u*]|.  The  signs  of  (Re  vT— t?  and  dm  y/l  — u*  are  found  by 
analytic  continuation  of  their  values  in  equation  (16),  as  the  cuts  of  Figure  2 
are  rotated  clockwise.  It  follows  that  sin  8 (Re  Vl— u5  + cos  8 dm  y/\  — u* 
is  positive  for  all  u on  the  contour  arg  u — —arg  k,  , and  consequent!;-  that 
eouation  (11)  approaches  zero  exponentially  as  | k,  j — »«>,  0 < arg  k,  < v. 
Moreover,  for  kx  real,  it  can  be  seen  that  ,(  — &,)  — &<(fci)),  the  bar  denoting 
complex  conjugate.  This  result  was  of  course  to  be  expected  since  k)  — 
&>o(k)  and  the  differential  operator  A*  *4*  k]  is  real  for  real  k,  . 

Using  Cauchy’s  theorem,  we  may  now  infer  rigorously  the  result  that  if 

a > 0,  6 > 0, 

(43)  4>.(-a  4*  ib)  = + *'&). 


We  write 
(44) 


<t>,(kt) 


J_  f dK<t>,{K) 
2ri  J K - kx  ’ 


the  contour  being  the  real  axis  from  -»  to  ® and  then  around  the  circle  at 
infinity  in  the  positive  kx  = K-plane.  Equation  (44)  is  justified  by  our  demon- 
stration that  <(>,  is  analytic  in  the  domain,  and  equation  (43)  follows  at  once  if 
we  not  ,*  that  the  integral  over  the  circle  at  infinity  vanishes  since  <f>,  has  been 
shown  to  become  zero  exponentially  on  this  circle. 

The  form  of  the  solution  so  far  obtained  is  awkward.  It  is  more  convenient 
to  introduce  2 J0  =*  Ho"  + H?'  as  in  Section  II.  This  is  not  possible  if  r =»  0 
since  Ho\0)  not  defined.  Restricting  ourselves  to  values'0  of  r > 0,  it  is  then 
easy  to  shew  using  Ho\kxru)  — — /fo‘)(A:,ru€,r)  that  equation  (17)  remains 
valid  and  is  an  analytic  function  of  k,  , provided  the  contour  C of  equation  (17) 
runs  from  infinity  to  the  origin  along  the  line  arg  u = r — arg  k , , and  then  from 
the  origin  to  infinity  along  the  line  arg  u = — arg  kx  . The  limits  within  which 
the  contour  can  be  deformed  can  be  inferred  from  equations  (21)  and  (22), 
which  remain  valid,  provided  8 is  interpreted  as  arg  kxu  instead  of  arg  u.  In 
III,  and  IV  then,  the  integrand  approaches  zero  exponentially  as  u — ♦<»  if 


“Obviously , in  justifying  equation  (17)  (cf.  Section  II)  we  should  have  stated  that  r > 0 
but  the  fact  that  our  formulas  involving  Hq1)  are  not  valid  for  r — 0 is  of  little  moment.  It  is 
important  however  that  a valid  solution  from  which  the  fields  at  r - 0 can  be  evaluated  for 
any  kx  is  available  (in  principle,  at  any  rate),  that  the  analysis  for  r > 0 yields  formulas  which 
can  be  understood  and  used  almost  everywhere,  and  that  there  is  no  reason  either  from  the 
physical  properties  or  from  the  appearance  of  the  formal  solution  to  expect  unusual  results  on 
the  line  r » 0. 
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(g  + h)  cos  $ + r sin  6 > 0,  or  it  may  be  deformed  from  its  original  arg  ktu  ■■  0 
anywhere  within  the  limits 

(45)  —tan-1  <argfr,+argu<T  — tan-1 * y 

In  I,  we  have  from  equation  (22) 

/>•  -i  * 4"  i v >■  ^ . . _j  t -t*  h 

(4tt)  tan  — - — < arg  t,  + arg  u < v 4-  tan  — - — . 

In  equations  (45)  and  (46) 

„ . , -i  z + h T 

0 < tan  < r. 

r 2 

The  inequalities  hold  up  to  and  including  the  limits;  i.e.,  < is  correct,  if  | fc,  \ > 0, 
as  can  be  seen  by  substituting  the  values  of  arg  k,v  at  the  limits  in  equations 
(21)  and  (22).  It  will  be  found  that  (Re  A (u)  is  not  aero  at  these  values  of  arg  klu, 
so  that  equation  (17)  converges  as  p — ♦ ® . 

These  results  will  be  used  in  the  following  section. 

VII.  Reception  of  a Pulse 

Suppose  the  wave  radiated  by  the  source  at  P (Figure  1)  is  given  by 

(47)  *»-g/(<-*) 

where  R is  the  distance  PS  and  /(z)  — 0 for  x < 0.  The  function  / (z)  satisfies 
sufficient  conditions  for  the  existence  of  a Laplace  transform,  but  is  otherwise 
arbitrary.  In  particular  /(z)  satisfies  the  condition  that  the  integral  from 
0 to  oo  of  | /(z)  | exists.  Then  /(z)  can  be  written  in  the  form 

(48)  J{x)  «=  ~ 

where  g(u)  is  a one- valued  regular  analytic  function  of  « in  the  entire  half  plane 
<m«>0. 

It  follows,  as  is  well  known,  that  the  solution  in  medium  1 is  given  by 
equation  (5)  where  now 

(49)  *,(t)  - ~ do,  e-<-‘g(„H>r(u) 

with$,(u)  « 4>r  of  eq.  (6).  The  functions $r(«)  have  been  derived  for  all  im  u > 0 
in  the  preceding  section.  In  particular  £,(«)  may  be  replaced  by  equations 
(9)— (11).  The  term  $,(u)  yields  a wave  which  is  received  at  a time  t — d/ci  and 
which  duplicates  in  shape  as  a function  of  time  the  shape  of  the  original  wave 
&,({).  This  is  also  a well  known  result,  and  is  demonstrated  by  the  usual  means 
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of  closing  the  contour  in  the  upper  w-plane.  We  still  have  to  examine  equation 

(60)  4>,(t)  « ^ J due~"“g(w)4>,(u) 

which  contains  all  the  deviations  from  geometrical  acoustics.  From  the  previous 
section  4>,  (u)  is  known  to  be  a regular  analytic  function  for  any  2,  r when  8m  « > 0. 
Since  g(u)  is  also  analytic  and  regular,  the  contour  may  be  deformed  at  will  in 
the  region  8m  u > 0. 

Deform  the  contour  so  that  it  makes  a small  half  circle  about  the  origin. 
This  is  done  to  avoid  the  difficulty  k,r  ■■  0.  Assume  r > 0.  Then,  in  equation 
(50),  $,(«)  is  given  by  eijuations  (17)  and  (43)  with  the  contour  C(kx)  defined  as 
in  the  previous  section,  lying  within  the  limits  of  equations  (45)  and  (46).  We 
shall  now  show  that  when  r < a the  integral,  equation  (50),  is  re ro  over  the 
infinite  half  circle  at  infinity  in  the  upper  w-plane,  8m  u > 0,  whenever 
t — d/ci  < 0.  Consequently  by  Cauchy’s  theorem  4>,(t)  — 0 whenever 
t — d/c,  < 0,  if  r < a. 

We  consider  first  the  quarter  circle  at  infinity  0 < arg  u < x/2.  For 
any  kx  in  this  region,  kx  = w/c,  , c,  real,  the  contours  C{kx)  can  be  deformed  to 
the  contour  C"  of  Figure  3.  In  order  that  this  deformation  be  possible,  according 
to  equations  (45)  and  (46),  it  is  only  necessary  that 

(61)  —tan-1  < arg  kx  + tan-1  T r < r — tan-1  r h 

' r ® z + h r 

and 

(52)  tan-*  - ^ < arg  £,+*■  — tan"*  — j— r < r + tan-1  - . 

Equations  (51)  and  (52)  are  satisfied11  provided  0 < arg  kx  < t/2. 

If  C,  represents  the  contour  along  the  quarter  circle  of  radius  p,  the  integral 
of  equation  (50),  integrated  over  this  contour,  can  now  be  written  as 


(53) 


<?(0  « Jc  duc~im,g(u)  du  exp  \ikx(h  + 2)  Vl  - u*| 


The  integral  over  C'  in  equation  (53)  converges  uniformly  as  u approaches 
infinity  for  any  uonC,.  C'  is  not  a function  of  w and  g(u)  is  bounded.  It  follows 
that  the  order  of  integration  can  be  interchanged  for  any  finite  p,  and  we  need 
only  show  that 


(54) 


exp  \ikx(h  -f  *)  vl  — ut\Hol\kxru)  « 0 


“They  tro  not  satisfied  if  arg  kx  > t/2;  this  is  the  reason  why  a line  of  reasoning  stemming 
from  equation  (43)  wi Vu  adopted. 
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As  p -+<»,  Ho''fJkxru)  can  be  replaced  by  its  asymptotic  expansion,  for 
any  u on  C.  uF(u)/VT - 1?  is  bounded.  The  limit  a a p — ♦»  is  therefore 
determined  by  the  behavior  of 

(55)  B(u)  = iAr,[A(u)  — c,f] 

with  A(u)  defined  by  equation  (19).  Letting  k,  - pe°,  it  is  evident  that  the 
limit  approaches  0 exponentially  as  p — ► «>  provided 


(56)  —sin  6 (Re  .l(u)  — cos  6 3m  vt(u)  + c,<  sin  0 < 0 


for  all  0 < 6 < w/2.  But  however,  in  Section  III  the  contour  C'  was  so  defined 
that  A(u)  — d + ip,  where  tj  > 0 for  aii  u on  C".  Consequently  equation  (56) 
is  satisfied1*  whenever  — d + cxt  <0.  On  the  other  hand  whenever  — d + clt>  0, 
there  will  be  a continuous  set  of  values  of  u on  C for  which,  for  a finite  range  of 
angles  6,  equation  (56)  is  not  satisfied.  It  follows,  using  equation  (43),  that  the 
contour  can  be  closed  in  the  upper  half  w-plane  whenever  t < d/cx  , and  cannot 
whenever  t>  d/cx . In  other  words  4>,{t)  of  equation  (49),  like  the  geometrically 
reflected  $,(f),  begins  at  t — d/cx  and  not  before. 

When  r > a the  contour  C(kx ) can  be  deformed  into  the  contours  C,  and 
C,  of  Figure  4 and  finally  into  the  contours  C[  and  C'2  of  Figure  5,  since  the  asymp- 
totic limits  of  C[  and  C2  are  the  same  as  C',  so  that  equations  (51)  and  (52) 
remain  relevant.  The  integrals  continue  to  be  uniformly  convergent,  and  the 
argument  goes  just  as  above.  The  integral  along  C2  leads  to  the  same  conclusion, 
namely  that  the  signal  is  received  when  t = d/cx  . The  integral  along  C[  leads 
to  equation  (56)  where  in  this  case  (Re  A(u)  - ra  + (z+h)  Vl  — a , so  that 
along  this  path  the  signal  is  received  only  when 


t 


ra  + (z  + h)  Vl  - q*  PA  + SB  AB 

Ci  C\  C) 


(Figure  1). 


That  no  poles  are  crossed  deforming  the  contours  when  r < a,  and  that  the 
pole  is  crossed  when  r > a only  when  a < 0 < 1 and  equation  (30)  is  satisfied 
is  shown  just  as  in  Sections  III  and  IV.  When  d>,(p)  contributes,  K0  is  replaced 
by  its  asymptotic  expansion  in  equation  (28)  and  it  is  found  that$,(p)  does  not 
contribute  until  cxt  equals  ( z + A)(l  - a*)ui/(l  - 0*)l/i.  That  this  time  lies 
between  the  times  of  appearance  of  the  head  wave  and  the  directly  reflected  wave 
follows  from  equation  (30).  Because  this  time  involves  0,  no  simple  geometrical 
interpretation  of  is  possible. 
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